Static and dynamical properties of elliptic quantum corrals 
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Due to their focalizing properties, elliptic quantum corrals present very interesting physical be- 
haviours. In this work we analyze static and dynamical properties of these systems. Results are pre- 
sented for realistic values of the parameters which might be useful for comparison with experiments. 
We study non-interacting corrals and their response to several kinds of external perturbations, ob- 
serving that, for realistic values of the Fermi level, the dynamics involves only a few number of 
excited states, making the system quite robust with respect to possible sources of decoherence. We 
also study the system in the presence of two S = 1/2 impurities located at its foci which interact 
with the electrons of the ellipse via a superexchange interaction J. This system is diagonalized 
numerically and properties such as the spin gap and spin-spin static and dynamical correlations are 
studied. We find that, for small J, both spins are locked in a singlet or triplet state, for even or odd 
filling respectively, and its spin dynamics consists mainly of a single peak above the spin gap. In 
this limit we can define an effective Heisenberg Hamiltonian to describe the low-energy properties. 
For larger J, more states are involved and the localized spins decorrelate, in a similar manner as 
the RKKY-Kondo transition for the two-impurity problem. 

PACS numbers: Pacs Numbers: 73.21.-b,75.75.+a,75.20.Hr 



INTRODUCTION 

The recent advancement in the knowledge and fabri- 
cation techniques of nanoscopic systems has revealed a 
great variety of novel physical behaviours, a striking ex- 
ample of which is the observation of the phenomenon 
now known as the quantum mirage This effect arises 
in nanoscopic systems called quantum corrals, where the 
electron wave functions are confined to a finite region by 
means of a barrier formed by atoms encircling a closed 
region in space. 

Quantum corrals are built by positioning atoms, typ- 
ically transition metal atoms, along a closed line on the 
clean surfaces of noble metals. In their recent experi- 
ments, Manoharan et. al.[l[ have built elliptical corrals 
with Co atoms on the (111) surface of Cu. The Cu (111) 
surface has a band of surface states, orthogonal to the 
bulk states, which can be represented as a two dimen- 
sional electron gas confined at the surface. The Fermi 
level is placed at 450 meV above the bottom of the sur- 
face state band. The atoms forming the corral act as 
scattering centers which tend to confine the surface elec- 
trons inside the corral. 

A number of theoretical papers analyzed these exper- 
iments with one or more impurities considering differ- 
ent configurations @, i i, [ESJzt Hi (see also 
complete review articles in [U 0|). In refs. [t| [1(3] 
it was suggested that, as a consequence of the focaliz- 
ing properties of quantum elliptic corrals, two impurities 
located at the foci of the system will strongly interact. 
Such a prediction has been supported by Stepanyuk et al. 
ll| who reported results of first principle calculations of 
the exchange coupling between magnetic impurities in- 



side quantum corrals. 

The physical properties of quantum corrals are affected 
by several factors, which include confinement, quantum 
interference and many-body effects. So far, interest has 
been focused on the static properties of quantum corrals, 
leaving the dynamical behaviour of such systems almost 
unexplored. 

In this work we have studied static and dynamical 
physical properties of elliptical quantum corrals for vari- 
ous models and perturbations, where the focalizing prop- 
erties of the system have important consequences in the 
physical behaviour. 

We have considered two cases. The first one consisted 
of an isolated ellipse describing a closed quantum sys- 
tem containing an arbitrary number of non-interacting 
electrons up to a certain Fermi energy. The system was 
perturbed by either local potentials located on the foci of 
the ellipse, or a potential in the form of a barrier cutting 
the ellipse through its minor axis. In the second case, 
two impurity spins interacting via a superexchange term 
with the electrons in the ellipse were added to the foci 
of the system. This many-body problem was treated nu- 
merically and static properties like the spin gap, and dy- 
namical responses, including real-time behaviour, where 
analyzed. 



NON-INTERACTING CASE 

Throughout this work we have adopted a simplified 
description of the system, based on a spinless tight- 
binding model defined on a discretized square lattice 
within the boundaries of the ellipse, with eccentricity 
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e = 0.6 and hopping element t* between nearest neigh- 
bours. The Hamiltonian in this case is simple and reads 
H e i = —t* Y^,Uj) c l c j' wnere the sites i and j run inside 
the ellipse which consists of around 1000 sites or atoms. 
Here we have considered a constant zero local energy 
in each site. As we will be analyzing low-energy wave 
functions with wave lengths much larger that the real 
interatomic distance, the underlying atomic distribution 
is irrelevant [Hj]. Our calculations could have also relied 
on the exact solutions for a continuum ellipse (Mathicu 
functions) . 

The Fermi energy was set to coincide with its 23 rd 
eigenstate IV'23), which is similar to the state at the Fermi 
energy in the experimental setup of Ref. [1] (see Fig. 
5(a)). This configuration is not crucial for the main con- 
clusions and other eccentricities and reference levels will 
have a similar behaviour, as long as the wave function 
bears an appreciable weight at the foci of the ellipse. 

We have studied the time evolution of the wave func- 
tion at the Fermi energy when, at time t = a constant 
negative or positive potential V\ is applied to focus 1: 
exp(— iHt)\i)2s), where H — H e i + Vi(t). Figs. 1, 2 and 
3 show the square of the amplitude of such a wave func- 
tion in foci 1 and 2 for t > 0, for three distinct cases 
corresponding to V\ = —1 (Fig. 1), V\ = 1 (Fig. 2) 
and V\ = V2 = — 1 (Fig. 3). We notice that in the two 
first figures, a finite time t elapses between the onset of 
the perturbation in focus 1 and the observation of the 
response in focus 2, as shown in the insets which is of 
the order of t ~ 20 H/2t* ~ 10~ 14 s, where we have used 
t* = 1 eV (see below). 

Another interesting point regards the oscillatory be- 
haviour of the two squared amplitudes for t > t$. The 
occurrence of oscillations in the probability amplitude in 
the presence of a time- independent potential (for t > 0) , 
as well as the overall "see-saw" or alternating behaviour 
between both foci, i.e., when the probability is high in 
one focus it is low in the other and viceversa (except 
for the higher order oscillations in focus 1 for negative 
V\), are intimately related to quantum interference ef- 
fects within the ellipse. The period of these oscillations 
is of the order of the energy difference corresponding to 
the main transitions involved, i.e. states 18, 23 and 29 
for this configuration. Using the parameters of Refs. 
0, 0, S] for the case of the ellipse on the copper [111] 
surface, i.e. distance between foci d — 70 A, total el- 
lipse area 85.7 nm 2 and t* = 1 eV, we obtain periods 
of the order of h/(E 2 3 — Ei$) — 5 x 10~ 14 s, where 
E23 — Eis = 0.0754 t* = 75 meV (the transition between 
states 23 and 29 has a similar value). We remark that 
this energy is one to two orders of magnitude larger than 
the Kondo temperature in quantum dots. An animated 
version of these behaviours is found in 14|. We notice 



potential, Fig. 2). 



that the amplitude of oscillation in focus 1 is larger than 
that in focus 2 when V\ is negative (attractive poten- 
tial, Fig. 1), and smaller when V\ is positive (repulsive 



0.16- 
0.12- 
[0.0B 
0.04 
0.00 



0.12-1 
0.Q8-I t 



50 100 150 




focus 2 



1000 1200 1400 1600 1800 2000 

tl(%I2t*) 

FIG. 1: Time dependence of the probability density in both 
foci for a negative potential, Vi = —1, applied in focus 1 at 
time t = 0. Insets depict short times and to is the time elapsed 
until the perturbation reaches the unperturbed focus. For 
t* = 1 eV, typical periods in this figure are around 5 x 10~ 14 s. 
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FIG. 2: Same as Fig. 1 but for a positive potential in focus 
1, Vi = l. 
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FIG. 3: Same as Fig. 1 but for a symmetric negative potential 
applied to foci 1 and 2, Vi = V2 = — 1. 

As a matter of corroboration, using the above set of 
parameters we obtain for the Fermi velocity of the elec- 
trons involved in the transmission of the information from 
one focus to the other the value vf = d/to = 0.8 x 10 8 
cm/s (were we have considered the focal distance d ~ 
70 x 10" 8 cto and t ~ 0.87 x 10~ 14 s), which is of the 
same order of magnitude as the one calculated in [H . 

We have also calculated the delay time and Fermi ve- 
locity for parameters corresponding to semiconductors, 
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where, for an effective mass m* = 0.07m e the energy dif- 
ference of the main transition is around leV and <o = 
0.06 x 10~ 14 s, which leads to a Fermi velocity of around 
10 9 cm/s and an electron density of 5.4 x 10 13 cm~ 2 , which 
is somewhat higher than the ty pica l carrier densities in 
semiconductor heterostructures 
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A particularly interesting case was the response to an 
infinitely large potential barrier along the minor axis ap- 
plied at times t > 0, starting from the same state as in the 
previous case (state 23). In Fig. 4(a) we see the evolution 
which is quite periodic with a frequency corresponding 
to the energy difference between the most relevant states 
involved (states 15, 16 and 25, 26 of the semi-ellipses, 
as shown in the approximate evolution where only these 
states were included). We, thus, observe here that due to 
the symmetry of the problem and the energy structure of 
the ellipse, only very few states take part and the time 
evolution becomes simple and robust. 
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FIG. 4: a) Time dependence of the probability density in focus 
1 for an infinite potential barrier applied along the minor 
axis at time t — 0. Solid and dot lines correspond to the 
time evolution of the exact and aproximate wave functions 
respectively, b) Modulus of the overlap between the non- 
perturbed and final wave functions after the potential barrier 
had been applied during an interval t\ and then removed (note 
the different parameters in the x-axis used both figures). 



INTERACTING CASE 

In addition to the one-particle behaviour studied in the 
previous section, we have analyzed the more interesting 
case of the effect of the focalizing properties of the ellipse 
on the interaction between two spins interacting via a 
superexchange term with the electron spins in the ellipse 
(Fig. 5). The Hamiltonian of the system reads: 



H = H el + J{S 1 -a 1 + S 2 -a 2 ) 1 



where 



Si ■ di = Sf 



-as: 



(i) 



(2) 



<t 4 + = c^c^, of = (riif — riii)/2, with the number 
operator and Cj CT the destruction operator of an electron 
with spin a in focus i of the ellipse. In the basis of eigen- 
states | a) of the ellipse, these local operators can be ex- 
panded as Cj CT = J2 a ^aiCaa where c a<7 and are the 
destruction operator and amplitude in state \a). H e i is 
the Hamiltonian of the isolated ellipse with infinite walls 
described in Sect. II. 

In this basis the spin operators are expressed as: 



<4 = Kri^C^Cc, 



2 1 



(3) 




b) 



> 



If, on the other hand, we leave the barrier potential 
on only for a finite time t\ and then calculate the over- 
lap of the original wave function to the remaining one 
after having removed the perturbation, we find another 
oscillatory behaviour between a nearly complete recov- 
ery of the original wave function and a nearly completely 
orthogonal state (see Fig. 4(b)). 



FIG. 5: a) 1 7/^23 1 2 for the non-interacting ellipse with eccentric- 
ity e = 0.6. White and black represent high and low electron 
densities respectively, b) Graphical representation of the sys- 
tem considered here where \a) represents an eigenstate of the 
non-perturbed ellipse. Si and S2 are the impurity spins at 
the foci of the ellipse which interact with the localized spins 
di via a superexchange interaction J. 
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Static properties 

We have solved this Hamiltonian numerically with ex- 
act diagonalization for a small number of levels and using 
the Lanczos technique for larger systems and different fill- 
ings corresponding to an even (closed shell) or odd (open 
shell) number of particles, N. The results for the spin 
gap A s are depicted in Fig. 6 where we find a differ- 
ent behaviour for even and odd number of electrons in 
the ellipse. For the former case we observe a slow in- 
crease of the spin gap as a function of the interaction 
parameter J which is quadratic and linear for small and 
large J respectively. For an odd number of electrons 
the behaviour is linear from the beginning and the gap is 
somewhat larger than for the even filling. One can under- 
stand this behaviour by using perturbation theory. For 
the even case (closed shell) the first correction to the en- 
ergy stems from second-order non-degenerate terms while 
for the odd filling (open shell), a first-order degenerate 
expansion is required, leading to a linear dependence. 
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FIG. 6: Behaviour of spin gap A s with J for even (triangules) 
and odd (circles) number of electrons in the ellipse. 



For a realistic value 
rameter for the Cu ellipse 



17] of the superexchange pa- 



J, 



= O.leV, we obtain 



J/2t* ~ 0.05 and a corresponding gap of around 10~ 5 
eV~ 0.12 K for even filling, temperature under which the 
system is in a stable singlet state. For the semiconduc- 
tor ellipse the spin gap could be around 10 times larger 
although we do not have any reliable estimation of J in 
this case. Thus, in both cases, copper or semiconducting 
corrals, this singlet state could be observed at sufficiently 
low temperatures. 

One of the most important consequences of impurity 
spins in the presence of an elliptical corral like the one 
considered here is their enhanced correlation when situ- 
ated at the foci of the system. For example, by analyz- 
ing the spin-spin correlation function in the ground state 
(V'oPi^lV'o) we find that they form a singlet or triplet 
for small J and even or odd particles, respectively (see 



Fig. 7(a)). A ferromagnetic interaction between local- 
ized spins for odd numbers of electrons in the ellipse was 
also found in [3] ■ For larger values of the interaction pa- 
rameter both impurity spins decorrelate. The results are 
weakly dependent on the number of levels considered, as 
also shown here. 
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FIG. 7: a) Spin-spin correlation between localized spins in 
the ground state for even and odd number of electrons and 
considering different number of levels in the ellipse, b) Top 
panel: same as (a), N even, for the first excited state (having 
total spin S — 1) showing the triplet character between the 
localized spins in both foci for small J. Lower panel: Corre- 
lation function between the localized spin Si and the spin of 
the itinerant electron o\ compared to the spin-spin correla- 
tion between localized spins for even TV (for the ground state 
in all cases). For small J the localized spins are entangled 
in a singlet state with no correlation with the electrons. For 
large J, these spins decouple while enhancing the local singlet 
interaction with the itinerant electron. 



It is also interesting to analyze the on-site correlation 
between the localized spin S\ and the spin of the elec- 
tron localized immediately above it, <j\. In Fig. 7(b) 
we show how this interaction, which is antiferromgnetic, 
strengthens by increasing J, while S\ decorrelates from 
Si. In addition, as J increases, the electrons in the ellipse 
become more localized near the impurity in order to take 
advantage of the magnetic interaction and the spin-spin 
correlation should grow towards the singlet value -0.75. 
We find a lower absolute value for large J and this is due 
to the finitcness of the system. In this limit, many more 
extended states are required to form the localized state 
(see Eq. 3). 
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Dynamical behaviour 

We have calculated the spin response function to a spin 
excitation performed in focus 1 [18J 



J=0.1 



where 



A(uj) = — — lim ImG(uj + irj + E ) 

7T Tj— >0+ 



G(z) = (^o|^(^- J ff)- 1 ^|^o) 



(4) 



(5) 
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FIG. 8: Spin spectral function ^4(w) for an even number of 
particles, three values of J and different number of levels in 
the ellipse showing a negligible "finite size" effect. For small 
values of J the spin gap is small and the first excitation is 
mainly a triplet state. See also Fig. 7(b), top panel. 
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FIG. 9: Same as Fig. 8 for odd N 
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FIG. 10: Time dependence of the spin-spin correlation func- 
tion between localized spins for the perturbed wave function 
IV'(i)} = U(t)S* \ipo) for small and moderate interactions. For 
small J the correlation is nearly a triplet as seen in previous 
figures. 



In Figs. 8 and 9 we show the results for different values 
of J and for even and odd fillings, respectively. For even 
N we find that the excitations for small J consist mainly 
of one peak located at the spin gap, which is small. In this 
case both localized spins form a very well defined singlet, 
whereas the first spin excitation has a triplet character 
as can be seen in Figs. 7(b) (top panel) and 10(a) (see 
below) . For larger values of J the spectral weight of this 
peak diminishes, other states get involved and the spec- 
tral function is more complicated (see also larger values 
of J in Fig. 7(b)). For odd ./V and small J both localized 
spins form a well defined triplet in the ground state and 
the first spin excitation corresponds to a spin flip of an 
electron spin, leaving S% and S2 still entangled in a triplet 
state. In both cases a negligible "size" dependence with 
the number of levels considered is found. 

The fact that, for small J and an even number of parti- 
cles, the first excitation is nearly purely a triplet can also 
be seen as a robust excitation which varies very slowly 
in time (see Fig. 10). Here we depict (ip(t)\SiS2\ip(t)} 
where \ip(t)) = U(t)Sf\ipo) and U(t) the time-evolution 
operator and \ipo) is the interacting ground state. For 
J = 0.1 we find a nearly constant behaviour, scrambling 
up with a larger amount of eigenstates for larger values 
of J. Note that the factor 1/4 which arises when com- 
paring this correlation to a pure triplet behaviour stems 
from the perturbation operator Sf applied to the ground 
state (the function \ip(t)) is not normalized). 
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CONCLUSIONS 



processes should be considered. 



We have studied static and dynamical properties of 
hard-wall elliptical quantum corrals subject to different 
kinds of perturbations. On one hand we considered non- 
interacting systems and applied localized perturbation 
potentials in the foci and along the minor axis of the el- 
lipse, finding robust oscillatory behaviours which are a 
consequence of the small amount of levels involved in the 
excitations. A finite elapse was found in the case that 
the potential was applied to one focus until the pertur- 
bation reaches the unperturbed one. On the other hand 
we studied the interacting case where two localized spins 
at the foci of the ellipse interacted antiferromagnetically 
with the itinerant electrons. In this case we found a dif- 
ferent behaviour for even and odd electrons which shows 
up in the spin gap, the character of the ground state and 
the excitations. For the even particle number case and 
small interaction J, both localized spins are entangled 
in a (quasi)-singlet state, forming a (quasi)-triplet in the 
first excitated state, which is the main spin excitation. 
These states are very robust, however, the spin gap is 
quadratically small with J. An interesting feature can 
be seen which resembles the RKKY-Kondo like transi- 
tion occurring in the two-impurity system when increas- 
ing the Kondo interaction J: While the localized spins 
form a singlet state for small interaction, this coupling 
decreases for larger J giving rise to an on-site singlet 
correlation between the spin and the itinerant electron. 

In this work we have shown that the focalizing prop- 
erties of elliptic corrals lead to an interesting behaviour 
between localized spins. For small J and even number 
of electrons in the ellipse, the system can be modelled 
by two spins interacting via an effective antifcrromag- 
netic interaction parameter H e ff — J e ffSi.S2, where 



J, 



eff 



A,. For odd number of electrons the interaction 



is ferromagnetic in the ground state. The character of 
this interaction can be controlled by changing the chem- 
ical potential of the system. 

For small to moderate values of J (J < 1) the 
main results synthesized above still hold in the case 
of more realistic models of quantum ellipses which 
include tunneling of the electrons in open corrals and 
inelastic processes with bulk electrons. In this parameter 
range the broadening of the relevant energy levels is 
smaller than their separation [l^, 13, [l9|, |2Q|. When 
larger interactions are included, higher levels which are 
more hybridized take part, and models including these 
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